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Transport coefficients associated with tire mass flux of impurities immersed in a moderately dense 
granular gas of hard disks or spheres described by the inelastic Enskog equation are obtained by 
means of the Chapman-Enskog expansion. The transport coefficients are determined as the solutions 
of a set of coupled linear integral equations recently derived for polydisperse granular mixtures [V. 
Garzo, J. W. Dufty and C. M. Hrenya, Phys. Rev. E 76, 031304 (2007)]. With the objective of 
obtaining theoretical expressions for the transport coefficients that are sufficiently accurate for highly 
inelastic collisions, we solve the above integral equations by using the second Sonine approximation. 
As a complementary route, we numerically solve by means of the direct simulation Monte Carlo 
method (DSMC) the inelastic Enskog equation to get the kinetic diffusion coefficient Do for two and 
three dimensions. We have observed in all our simulations that the disagreement, for arbitrarily 
large inelasticity, in the values of both solutions (DSMC and second Sonine approximation) is 
less than 4%. Moreover, we show that the second Sonine approximation to Do yields a dramatic 
improvement (up to 50%) over the first Sonine approximation for impurity particles lighter than 
the surrounding gas and in the range of large inelasticity. The results reported in this paper are 
of direct application in important problems in granular flows, such as segregation driven by gravity 
and a thermal gradient. We analyze here the segregation criteria that result from our theoretical 
expressions of the transport coefficients. 



I. INTRODUCTION 



The theoretical basis for a hydrodynamic description of ordinary (elastic) gases is well established at low density 
using the Boltzmann kinetic equation. However, for moderately dense gases there is no accurate and practical 
generalization of the Boltzmann equation except for the idealized hard sphere fluid. For this intermolecular potential, 
the Enskog kinetic equation takes into account the dominant positional corrections to the Boltzmann equation due to 
excluded volume effects but, like the Boltzmann equation, neglects velocity correlations (molecular chaos assumption) 
among particles which are about to collide The Enskog equation is the only available theory at present for 
making explicit calculations of the transport properties of moderately dense gases and in any case the molecular chaos 
assumption is expected to fail only in much denser systems (solid volume fractions typically are larger than 0.4) d]. 
The extension to mixtures requires a revision of the original Enskog theory for thermodynamic consistency (revised 
Enskog theory, or RET) [|, and its application to hydrodynamics and Navier-Stokes (NS) transport coefficients was 
carried out more than twenty years ago 

Early attempts [1, d, 0, [1] to extend the study of Lopez de Haro, Cohen and Kincaid (3] to inelastic hard sphere 
mixtures were restricted to nearly elastic systems. In this case, the effect of inelasticity in the collisions is taken into 
account only by the presence of a sink term in the energy balance equation and. as a consequence, the expressions 
of the NS transport coefficients arc the same as those obtained for clastic systems. Moreover, those early works also 
assume energy equipartition and so the partial temperatures for each species are equal to the granular temperature. 
However, as the dissipation increases, different species of a granular mixture have different partial temperatures 
Ti and consequently, the energy equipartition is seriously broken down (Ti ^ T) [^, [lol . [ll| . The failure of energy 
equipartition in granular fluids has also been confirmed by computer simulations [12| and observed in real experiments 
[l3| of agitated mixtures. Results in the literature show that the deviation from equipartition depends on the size 
and mass ratios of the particles of each species and the coefficients of restitution of the system. 

A mo re g eneral extension of the RET to inelastic collisions has been recently carried out by Garzo, Dufty and 
Hrenya 14, [l^. This theory covers some of the aspects not taken into account in previous works f^, @, 0, [g and 



extends previous results derived for monodisperse dense systems [l6l.[T7j and dilute binary mixtures [l8| |. Specifically, 
(i) it goes beyond the weak dissipation limit so that it is expected to be applicable to a wide range of coefficients 
of restitution, (ii) it takes into account the non-cquipartition of granular energy, and (iii) it has been formulated 
for multicomponcnt systems without limits on the number of components. Therefore, this theory [3 subsumes all 
previous analysis for both ordinary and granular gases, which are recovered in the appropriate limits ^ [l6l. [T^. [l8| . 
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Nevertheless, as in the elastic case [J], although the results are exact in the first order of spatial gradients, the explicit 
form of the NS transport coefficients requires to solve a set of linear integral equations. The standard method to 
get the kinetic and collisional contributions to transport coefficients and cooling rate consists of approximating the 
solutions to these integral equations by Maxwellians (at different temperatures) times truncated Sonine polynomial 
expansions. For simplicity, usually only the lowest Sonine polynomial (first Sonine approximation) is retained [15| . and 
the results obtained from this approximation compare very well with Monte Carlo simulations of the Enskog equation 
in the case of the shear viscosity coefficient of a mixture heated by an external thermostat [T9| . However, exceptions 
to this good agreement are extreme mass or size ratios and strong dissipation, although these discrepancies could be 
mitigated in part if one considers higher-order terms in the Sonine polynomial expansion, as previous studies in the 
dilute limit indicate [20| . In fact, recent works for monodisperse gases have shown that higher order terms in Sonine 
polynomial expansions become increasingly important in the range of moderate and strong inelasticities [2^ . [23l . and 
for this reason it has been of interest to calculate transport coefficients with more refined Sonine approaches [2J, [2^. 
However, the above works have been mainly focused with low density granular gases and many of the problems of 
practical interest in granular gases lie in the range of moderate densities. For this reason, it is important to determine 
the degree of accuracy of at least the first Sonine approximation compared to the second Sonine approximation for 
dense granular gases. Therefore, by testing an eventual gain of accuracy with higher order Sonine approximations, 
our results in the present work will contribute to the debate in the literature on the validity of a hydrodynamic 
description of granular gases [2l[. In addition, the range of high inelasticities has growing interest in experimental 
works [2^, [l^l . Another motivation to improve the evaluation of the NS transport coefficients lies in the fact that 
the reference homogeneous cooling state (HCS) is known to suffer a clustering instability, with inter-cluster distance 
inversely proportional to inelasticity psj . In this context, we believe that a more accurate description of the HCS in 
the range of mild and strong inelasticities may help to refine the understanding of this interesting instability. 

Needless to say, the evaluation of the NS transport coefficients for a dense granular mixture beyond the first 
Sonine approximation is quite intricate, due mainly to the coupling among the different integral equations obeying 
the transport coefficients. We will thus make a first approach to the problem by considering the simple situation of 
a granular binary mixture where the concentration of one of the species (of mass mg and diameter (Tq) is very small 
(impurity or tracer limit). Moreover, the tracer limit has been of much interest in recent literature, for example in 
granular segregation problems [2^ [30, HH, H^l ■ In the case of a tracer immersed in a dense granular gas, and as in a 
previous study for dilute gases [20|, one can assume that (i) the state of the dense gas (excess component of mass m 
and diameter a) is not affected by the presence of impurities or tracer particles and, (ii) one can also neglect collisions 
among tracer particles in their corresponding kinetic equation. As a consequence, the velocity distribution function 
/ of the gas verifies a closed Enskog equation while the velocity distribution function /o of the tracer particles obeys 
the linear Enskog-Lorentz equation, which greatly simplifies the development of Chapman-Enskog theory. 

Under the above conditions, since the pressure tensor and heat fiux of the mixture (gas plus impurities) is the same 
as that for the gas [H, [l3| , the mass transport of impurities jg is the relevant flux of the tracer problem. To first 
order in the spatial gradients, three transport coefficients are involved in the constitutive equation for jg: the kinetic 
diffusion coefficient Dp, the mutual diffusion coefficient D and the thermal diffusion coefficient . Thus, the mass 
flux jg has the form 

jg = ^i?oVno °-DVn - ^D^VT, 1 

P P T 

where p — mn is the total mass density, rig is the number density of the impurities, n is the number density of the gas 
particles, and T is the granular temperature. Therefore, the main goal of this paper is to determine Dq, D and I?^ up 
to the second Sonine approximation in terms of the coefficients of restitution for the impurity-gas (ag) and gas-gas 
(a) collisions, the parameters of the system (masses and sizes) and the solid volume fraction </> occupied by the gas. 
The calculations are rather intricate and we have taken advantage of some previous calculations performed in Ref. [l5| 
for multicomponent systems. In particular, a previous expression for the coefficient Dq obtained in the second Sonine 
approximation for a dilute gas {(f) ~ 0) is recovered [l^l • Analogously to the previous analysis of the shear viscosity 
coefficient [l^ , kinetic theory predictions for the diffusion coefficient Dq are compared with numerical solutions of the 
Enskog equation by using the well-known direct simulation Monte Carlo (DSMC) method [11]. In the simulations, 
the diffusion coefficient is computed from the mean-square displacement of impurities immersed in a dense granular 
gas undergoing the homogeneous cooling state [20| . Although the problem is time-dependent, a transformation to a 
convenient set of dimensionless time and space variables 34] allows one to get a stationary diffusion equation where 
the coefficient Dq can be measured in the hydrodynamic regime (times large compared with the characteristic mean 
free time). 

Finally, once the explicit expression of the transport coefficients associated with the mass flux arc at hand, a 
segregation criterion based on thermal diffusion is derived. This criterion shows the transition between the well- 
known Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE) by varying the different parameters of the 
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system. This study complements a previous analysis recently carried out by one of the authors [32j for a driven 
(heated) dense gas. As expected, our results show that the form of the phase-diagrams for the BNE/RBNE transition 
depends sensitively on the value of gravity relative to the thermal gradient and so it is possible to switch between 
both states for given values of the parameters of the system. 

The plan of the paper is as follows. In Sec. [TTl we describe the problem we are interested in and offer a short 
summary of the set of inelastic Enskog equations for the gas and the impurities. Section Hill deals with the application 
of the Chapman-Enskog method [sH] to solve the Enskog-Lorcntz equation and get the set of coupled linear integral 
equations verifying the transport coefficients Dq, D and . Then, these integral equations are approximately solved 
up to the second Sonine approximation. Some technical details of the calculations are given in Appendices [X] and |B] 
In Sec. llVl we illustrate the dependence of the transport coefficients on the parameters of the system and compare the 
theoretical results for the coefficient Dq obtained in the first and second Sonine approximation with those obtained by 
means of Monte Carlo simulations of the Enskog-Lorentz equation for disks {d = 2) and spheres (d = 3). Segregation 
by thermal diffusion is studied in Sec. |V] and the paper is closed in Sec. IVII with a brief discussion on the results 
derived. 



II. DESCRIPTION OF THE PROBLEM 



Let us consider a binary mixture of inelastic particles with collision rules according to the smooth hard sphere 
model. Our system is described by the revised Enskog kinetic equation [36l. [37j and, as a starting point, we consider 
in this work the special case where the concentration of one of the components (the tracer or intruder) is very small 
compared to that of the other (solvent or excess) component. In this limit, the state of the granular gas (the solvent) 
is not affected by the presence of the tracer particles and also the mutual interactions among the tracer particles 
can be neglected as compared with their interactions with the particles of the solvent. At a kinetic theory level, this 
implies that the velocity distribution function of the solvent (/) and of the tracer particles (/o) obey respectively the 
closed (nonlinear) Enskog equation and the (linear) Enskog-Lorentz equation. This is formally equivalent to study an 
impurity or intruder in a dense granular gas, and this will be the terminology used here. Since in the tracer limit the 
pressure tensor and the heat flux of dense binary mixtures are the same as those of the pure excess component (in 
the absence of the tracer), here we will be interested in the evaluation of the transport coefficients defining the mass 
flux of the intruder. 

Let us start by offering a short review on some basic aspects of the set of inelastic Enskog equations for the gas 
and the intruder. The granular dense gas is composed by smooth inelastic hard disks {d = 2) or spheres (d = 3) of 
mass m and diameter a. The inelasticity of collisions among all pairs is accounted for by a constant coefficient of 
normal restitution a (0 < a < 1) that only affects the translational degrees of freedom of grains. The granular gas 
is in the presence of a gravitational field g = —gSz, where g is a positive constant and is the unit vector in the 
positive direction of the z axis. At moderate densities, we assume that the time evolution of the one-particle velocity 
distribution function of the gas /(r,v, <) is given by the Enskog equation (36l. 

5t + V • V + g /(r, V, t) = J[v|/(i), f{t)l (2) 

where the Enskog collision operator J[v|/, /] is 

J[ri,vi I /(i),/(t)] = a^-' Jd^,2 I d?e(a-gi2)(a-gi2) 

X [a-^x{ri,ri - (t) f{ri,VI;t)f{ri - (T,V^;t) 

-X ri +cr) /(ri,vi;i)/(ri +cr,v2;i)] ■ (3) 

Here, a is the hard sphere diameter, S is a unit vector along their line of centers, O is the Hcaviside step function, 
and gi2 = Vl — V2 is the relative velocity. The primes on the velocities denote the initial values {v'/ , V2 } that lead to 
{vi,V2} following a binary collision in the hard sphere model: 

v" = Vl - i (f + a-^) {a ■ gi2)?, v^' = V2 + J (1 + a-^) {a ■ gi2)?. (4) 

The quantity x(ri,ri +(t | n{t)) is the pair correlation function at contact as a functional of the nonequilibrium 
density field n{r,t), where 

n(r,t) = J dv/(r,v,0. (5) 
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In addition, the flow velocity and the granular temperature arc defined respectively as 

u{r,t} = / ,N / rfvv/(r,v,i), (6) 
n(r,t) J 



(7) 



where V(r,i:) = v — \i[r,t) is the peculiar velocity. The macroscopic balance equations for number density n, 
momentum density mu, and energy density ^nT follow directly from Eq. ^ by multiplying with 1, mv, and ^mv^ 
and integrating over v: 

Dtn + nV • u , (8) 
DtU+{mn)-^V -P = g, (9) 

AT + (V • q + P., V,u,) = -CT , (10) 

where Dt = 9t + u • V is the material time derivative. The microscopic expressions for the pressure tensor P, the 
heat flux q, and the cooling rate ( in terms of the velocity distribution function / can be found in Refs. [lB| and 
Their explicit forms will be omitted here for brevity. 

Let us suppose now that an impurity or intruder of mass mp and diameter ctq is added to the gas. As said before, 
the presence of the intruder does not have any effect on the state of the gas, so that its velocity distribution function 
is still determined by the Enskog equation In addition, the macroscopic flow velocity and temperature for the 
mixture composed by the dense gas plus the intruder are the same as those for the gas, namely they are given by Eqs. 
^ and d?]), respectively. Under these conditions, the velocity distribution function /o(r, v, i) of the intruder satisfies 
the linear Enskog-Lorentz equation 

+ V • V + g • 7^ ./o(r, V, t) = JoMMt)Jit)], (11) 



dt 9v 
where the collision operator Jo[v|/o(t), f{t)] is now 

Jo[ri,vi|/o(t),/(t)] = a""-' J dw2 J d?e(? • gi2)(? • gia) 

X [ao'^Xo (ri,ri - ff)/o(ri,v'/;t)/(ri - o^,v^';t) 

-Xo (ri,ri + o')/o(ri,vi;t)/(ri + W,V2;t)] . (12) 

Here, W = Wa, W = (ctq + cr)/2, ao (0 < ao < 1) is the coefficient of restitution for intruder-gas collisions, and xo is 
the pair correlation function for intruder-gas pairs at contact. The precollisional velocities are given by 

v" = vi ^ (l + a^^) {a ■ gi2)S, 

mo + m 

v'i = V2 + (1 + ao"') • gi2)^. (13) 

Too -I- to 



As shown in Ref. [llj, the operator Jq [v|/o, /] is the same as that of an elastic impurity (ao = 1) with an effective 
mass 

'niQ=mo + - ^(toq-I-to). (14) 

1 + ao 

The number density for the intruder is 

no(r,t) = / (iv/o(r,v,t). (15) 
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The intruder may freely loose or gain momentum and energy in its interactions with the particles of the gas and, 
therefore, these are not invariants of the collision operator Jo[v\fo, f]- Only the number density uq is conserved, 
whose continuity equation is directly obtained from Eq. (fTTj) 

Ano+noV-u+^^ = 0, (16) 
mo 

where jo is the mass flux for the intruder, relative to the local flow u, 

jo = mo J dvV fo{r,v,t). (17) 
At a kinetic level, an interesting quantity is the local temperature of the intruder, defined as 

^°(^'^)-^/ ^-^^/o(r,v,t). (18) 

This quantity measures the mean kinetic energy of the intruder. As confirmed by computer simulations [l3|, experi- 
ments [Tsl and kinetic theory calculations Q, the global temperature T and the temperature of the intruder To a-re 
in general different, so that the granular energy per particle is not equally distributed between both components of 
the system. 

III. MASS TRANSPORT OF IMPURITIES 

In order to compute the mass flux jo of impurities to first order in the spatial gradients, we solve the Enskog-Lorentz 
equation by means of the Chapman-Enskog (CE) expansion 1351 . This method, nowadays extensively used and tested 
in a variety of problems in the field of rapid granular flows [S^ , assumes the existence of a normal solution in which 
all the space and time dependence of fo occurs through the hydrodynamic fields no, n, u and T. The CE procedure 
generates the normal solution explicitly by means of an expansion in the gradients of the fields: 

/o=/r+e/(i''+--- , (19) 

where e is a formal parameter measuring the nonuniformity of the system. The application of the CE method to the 
Enskog equation for polydispersc granular mixtures has been carried out very recently p^. Il5|. Using those results, 
we consider here the tracer limit {xo = no/n — > 0) for the linear integral equations defining the transport coefficients 
associated with the mass flux. The flrst-order contribution jg^^ to the mass flux is given by Eq. ([1]), where the kinetic 
diffusion coefficient Do, the mutual diffusion coefficient D, and the thermal diffusion coefficient are defined as 

D^--^ JdvV-A,{V), (20) 
JdvV-Bo{V), (21) 



7710 710 d 



D 



iydvV.Co(V). (22) 



The quantities Aq (V) , Bo (V) , and Co (V) are the solutions of the following set of coupled linear integral equations 

ic(°)^ ■ (VA,) - Ic^'^Ao JnA„/(°)] = Ao + 4"\C,A], (23) 
^C(°)^-(VBo)-Jr[Bo,/(")]=Bo, (24) 
ic(°)^ ■ (VCo) - n^A, - Jf [Co, = Co + Jr ,C], (25) 
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where (■'■"^ is the coohng rate to zeroth order (local homogeneous cooling state) and Jq ' [X, Y] is the operator 
[v,\X,Y] = x^^'"' / J dae{a ■ gi2)(S • gi2) [a^^Xmrm - X(Vi)y(V2)] , 



(26) 



where Xo"'' i^ the intruder-gas pair correlation function at zeroth order. The inhomogeneous terms of the integral 
equations ([23l) -([25 |) are defined by 



A.fv/("; 



(27) 



Bo (V) = -V/^ 



(0) 



(28) 



T,no 



(0) 



l + 2'^-V"V(l + a) 



In Eqs. ([27)) -((29 |) . the pressure p is given by 

p = nT 

uj = ao/a is the size ratio and /xo is the chemical potential of the intruder. Furthermore, 



^ 2^-idr((i/2) 

is the solid volume fraction and the operator /Co,i[^] is defined as 

iCoAx]^a^x^^^ I dv2 / d5e(^ •gi2)(S-gi2)i?. kvr(v'/)x(v^') + /r(Vi)x(V2) 



(29) 



(30) 



(31) 



(32) 



Upon writing Eqs. use has been made of the expression of the first-order distribution function f^^^ of 

jas particles. Its form has been derived in Refs. [IB, and reads 



/(I) = A • VT + C • Vn + X> : Vu + • u. 



(33) 



where the coefficients A.. CD and E are functions of the peculiar velocity V and the hydrodynamic fields. According 
to Eqs. ([23|) - ((25)) . only the coefficients A. and C are involved in the evaluation of the transport coefficients Z?o, 
D and . The linear integral equations verifying A and C as well as their expressions up to the second Sonine 
approximation are given in Appendix \K\ 

It is worthwhile to remark that the quantities .4fl and Co associated with the intruder are coupled with their 
corresponding counterparts A. and C of the host gas through the integral equations ((23)) and (|25|) . respectively. A 
direct consequence of this coupling is that the mass flux of the intruder ([T]) inherits gradient terms (Vn and VT) from 
those of the autonomous host equations. Moreover, the external field does not occur in the constitutive equation ([T]) 
for the mass flux. This is due to the particular form of the gravitational force. 



IV. SECOND SONINE POLYNOMIAL APPROXIMATION 



For practical purposes, the integral equations (|23[) - (|25p can be solved by using a Sonine polynomial expansion. 
With the motivations explained in the Introduction, our goal here is to determine the diffusion coefficients Dq and D 
and the thermal diffusion coefficient I?^ up to the second Sonine approximation. In this case, the quantities Aq, 3q 
and Co are approximated by 



A(V) ^ -.fo,M(V) 



YD^ + aoSoiV) 



(34) 



5o(V)^-/o,m(V) 



-^VDo + boSoiY) 

/OJO 



(35) 
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Co(V) ^ -/o,m(V) 

where 



^YD + coSoiV) 

noTo 



(36) 



So(V) = ( ImoV^ - ) V, (37) 



2 2 

and /o,j\f(V) is a Maxwcllian distribution at the temperature To of the intruder, i.e., 



/o^/W^nofeV'-pf-^). (38) 



The coefficients oq, 6o and cg are defined as 



«o = -377^^ / '^v So(V) • A)(V), (39) 



d{d + 


2) noTi 


2 


mo 


d{d + 


2) noT^ 


2 


mo 



y dv So(V) • So(V), (40) 



Co = / rfvSo(V)-Co(V). (4f) 



The transport coefficients -Doi D and -D"^ as well as the second Sonine coefficients Cq, &o ^-nd Co are determined by 
substitution of Eqs. into the integral equations ([^ - (PS)) . multiphcation of these equations by moV and by 

So(V), and integration over velocity. The details are carried out in Appendices IB] and [Cl and only the final expressions 
will be presented here. 

The second Sonine approximations -Do [2], D[2] and -D^[2] can be written, respectively, as 

Do[2] =^^(a,ao,mo/m,(7o/(T,<^)i?o[l], (42) 
D[2] = G(a, ao, mo/m, ao/a, 0)Z?[f ], (43) 

D^[2] ^ H{a,ao,mo/m,ao/a,^)D^[l], (44) 

where F, G and H are nonlinear functions of the mass and size ratios, the coefficients of restitution and the solid 
volume fraction. The explicit forms of G and H are given by Eqs. (jBf 9p . (|B24p and (jBf 7[) . respectively. In Eqs. 
(|4T|) - ([44|) . Do[l], D[l\ and -D"^[f] refer to the first Sonine approximations to Do, F> and D"^, respectively. Their explicit 
expressions were already determined in Ref. [ist for arbitrary composition. In terms of the transport coefficients, the 
new calculations in the present work are the functions F, G, and H. In the tracer limit {xq 0), the expressions of 
Do[l]; D[l] and -D"^[f] reduce, respectively, to 

Do[l] = (45) 

m = ^-^, (46) 
moJ^ i^l - 

^"[1]-^^. (47) 



Here, ly ~ na^^ ^ ^2T jm is an effective collision frequency, 7 ^T^jT is the temperature ratio, 

/-(O) /2_(d-l)/2 

C = — = Y(°)(I-a2) (48) 

^ V dT{dl2) ^ " ^ 
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FIG. 1: (color online) Plot of the reduced kinetic diffusion coefficient Do{a) / Do{l) as a function of the (common) coefficient 
of restitution a = ao for the systems (mo/m = 4, ao/cr — 2) (a) and (mo/m = 0.5, cro/cr = 0.8) (b) in the case of a three- 
dimensional gas with (f> — 0.1. The solid lines correspond to the second Sonine approximation while the dashed lines refer to 
the first Sonine approximation. Here, -Do(l) is the elastic value of the kinetic diffusion coefficient consistently obtained in each 
approximation. 
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FIG. 2: (color online) Plot of the reduced mutual diffusion coefficient D{a)/ D{1) as a function of the (common) coefficient 
of restitution a = ao for the systems (mo/m = 4, ao/cr — 2) (a) and (mo/m = 0.5, cro/cr = 0.8) (b) in the case of a three- 
dimensional gas with (f) = O.f. The solid lines correspond to the second Sonine approximation while the dashed lines refer to 
the first Sonine approximation. Here, D{1) is the elastic value of the kinetic diffusion coefficient consistently obtained in each 
approximation. 
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FIG. 3: (color online) Plot of the reduced thermal diffusion coefficient {a)/D^(l) as a function of the (common) coefficient 
of restitution a = oq for the systems (mo/m — 4, cro/cr — 2) (a) and (mo/m = 0.5, cro/cr = 0.8) (b) in the case of a three- 
dimensional gas with (j) = O.f. The solid lines correspond to the second Sonine approximation while the dashed lines refer to 
the first Sonine approximation. Here, (1) is the elastic value of the thermal diffusion coefficient consistently obtained in 
each approximation. 
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FIG. 4: (color online) Plot of the reduced kinetic diffusion coefficient Do{a) / Do{l) as a function of the (common) coefRcient 
of restitution a = ao for mo/m = 1/8, ao/a — 1/2 and (j) — 0. The left panel is for hard spheres (d = 3) while the right 
panel is for hard disks (d = 2). The solid lines correspond to the second Sonine approximation, the dashed lines refer to the 
first Sonine approximation and the dotted lines are the modified Sonine approximation. The symbols are the results obtained 
from Monte Carlo simulations. Here, -Do(l) is the elastic value of the thermal diffusion coefRcient consistently obtained in each 
approximation. 
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FIG. 5: (color online) Plot of the reduced kinetic diffusion coefficient Do{a) /Do{l) as a function of the (common) coefficient 
of restitution a = ao for mo/m = 1/5, ao/a — 1/2 and = 0.2. The left panel is for hard spheres (d = 3) while the right 
panel is for hard disks {d — 2). The solid lines correspond to the second Sonine approximation, the deished lines refer to the 
first Sonine approximation and the dotted lines are the modified Sonine approximation. The symbols are the results obtained 
from Monte Carlo simulations. Here, -Do(l) is the elastic value of the thermal diffusion coefficient consistently obtained in each 
approximation. 



is the (reduced) cooling rate and the reduced quantities X^, Yj' and i^l are given by Eqs. (jBlSp . (|B2ip and (jCip . 
respectively. 

In general, the first and second Sonine approximations for the transport coefficients of the mass flux have a complex 
dependence on the coefficients of restitution, the solid fraction and the mass and size ratios. Thus, before analyzing 
this dependence, it is instructive to consider some special limits. In the elastic limit (a = an = 1) of a three- 
dimensional system, one recovers previous results for a gas mixture of elastic hard spheres [Sa . |40| . Moreover, 
in the case of mechanically equivalent particles (toq = rn, (Jq = a, a = ao), as expected, one gets D"^[2] — 0, 
Do[2] ~ —{m/xomo)D[2], and so 

jW ^ ^!^Dom\/xo, (49) 

where xq = no/n is the mole fraction of impurities. Moreover, in the case of a dilute gas ((/) = 0), the expression of the 
kinetic diffusion coefficient Dq [2] coincides with the one previously derived by one of the authors (20| by assuming that 
the solvent is in the homogeneous cooling state. All these results confirm the self-consistency of the results reported 
in this paper. 
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FIG. 6: (color online) Plot of the reduced kinetic diffusion coefficient Do{a) / Do{l) as a function of the (common) coefficient 
of restitution a = ao for mo/m — 2, ao/cr — 2 and = 0.2. The left panel is for hard spheres {d = 3) while the right panel is 
for hard disks (d = 2). The solid lines correspond to the second Sonine approximation and the dashed lines refer to the first 
Sonine approximation. The symbols are the results obtained from Monte Carlo simulations. Here, -Do(l) is the elastic value of 
the thermal diffusion coefficient consistently obtained in each approximation. 



SOME ILLUSTRATIVE EXAMPLES AND COMPARISON WITH MONTE CARLO SIMULATIONS 



The expressions for the (reduced) transport coefficients Do[2]/Do[l], D[2]/D[l] and D'^ [2] / D'^ [1] depend 

on many parameters: {mo/m, cto/ct, a, ao, </)}, or equivalently [ll| {mo/m, ctq/ct, a, 0}, where mj is defined by Eq. 



p4)) . This complexity exists in the elastic limit as well [J], except for the dependence on the coefficients of restitution. 
Thus, to show more clearly the influence of dissipation on the transport coefficients we normalize them with respect 
to their values in the elastic limit. Also, for simplicity, we take the simplest case of common coefficient of restitution 
a = ao so that, the parameter space has four independent quantities: {mo/m, cto/ct, a, 4>}. 

In order to get the explicit dependence of the transport coefficients on the above four parameters one has to give 
the forms of the pair correlation functions x^^'' Xtl' three-dimensional case (d = 3), a good approximation 

for x^"-* is provided by the Carnahan-Starling form [4l| 



X 



(0) 



1 



(1- 



while the intruder-gas pair correlation function is given by [42 

1 



(0) 



1 



3 



(1+C^)2 (1 



(50) 



(51) 



where we recall that uj = (To/cr is the diameter ratio. The expression for the chemical potential of the intruder 
consistent with the approximation ([5T|) is [i^ 



^ = C3+lnno-ln(l-0)-f 3c^- 



ln(l 



0(2 -</>)' 



(1 



21n(l - 



(1-60 + 3(^2) 



(1-0)3 

where C3 is a constant. For a two-dimensional gas (d = 2), x^^^ ^ind Xq"' are approximately given by [44 

AO) 



1 



16^ 



(1-0)2 



(52) 



(53) 



(0) 

Xo 



1 



The chemical potential is now given by [45 



— = C2 + In no - ln(l - 



— UJ 

4 



90 



1 +LJ (1 - 



ln(l - 



—LlJ 

8 



(54) 



(55) 
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where C2 is a constant. 

In Figs. [T]-[31 we plot the transport coefficients for inelastic hard spheres {d = 3) as functions of the coefficient of 
restitution for two different systems. Each transport coefficient has been reduced with respect to its elastic value 
consistently obtained in each approximation. The dashed lines refer to the first Sonine approximation while the 
solid lines correspond to the second Sonine approximation. We observe that in general the first Sonine polynomial 
approximation quantitatively differs from the second Sonine approach as the dissipation increases for sufficiently 
small values of the mass ratio mo/m and/or the size ratio at^ja. For these cases, the first Sonine approximation is 
not sufficient to capture the influence of dissipation on mass transport. However, the predictions of the first Sonine 
correction improve significantly as mo/m and/or cro/fj increases so that the former accurately describes the mass 
transport in this range of values of the mass and size ratios, even for strong inelasticity. These findings on the 
convergence of the Sonine polynomial expansion are quite similar to those obtained for elastic systems [40[ and for 
granular gases at low-density [20| . 

To check the reliability of the first and second Sonine approximations, we have performed Monte Carlo simulations 
of the Enskog equation for the dense granular gas with the tracer particles. We have extracted from these simulations 
the kinetic diffusion coefficient Dq of impurities in a granular dense gas in the homogeneous cooling state (HCS). 
This coefficient can be obtained from the mean square displacement of the intruder particle after a time interval t as 

|(|r(t)-r(0)P) = ^, (56) 

where |r(t) — r(0)| is the distance travelled by the intruder from t = until time t. Equation ([56|) is the Einstein form 
of the diffusion coefficient. This relation can be used also in Monte Carlo simulations of granular gases to measure the 
diffusion coefficient (see for example a previous work on transport of impurities in a dilute granular gas in Ref. [20|). 
In an unbounded system like ours, the DSMC method has two steps that are repeated in each time iteration: the first 
step takes care of the particles drift and the second step accounts for the collisions among particles. The extension 
of the DSMC method to study the diffusion of impurities in a dense granular gas in the HCS requires the changes 
~^ xJ[fif] ^-nd Jo[fo, f] XoJoifo, f]- For the DSMC method to work appropriately, the time step needs 
to be small in comparison with the microscopic time scale of the problem (which is set by the collision frequency ly) 
and wc also need a sufficiently high number of simulated particles [s^. Thus, we have used in the simulations of this 
work a time step St = 2.5 x 10~^i/~^ and N = 2 x 10^ simulated particles for each species [i^. To our knowledge 
and since we are interested in the complete range of values of a, we present the first DSMC data on dense granular 
gases for coefficients of restitution as low as a = 0.1. More details on the application of the DSMC method [3^ to 
this diffusion problem can be found in Ref. [l^l . 

If a hydrodynamic description (or normal solution in the context of the CE method) applies, then the diffusion 
coefficient Do(t) depends on time only through its dependence on the temperature T{t). In this case, after a transient 
regime, the reduced diffusion coefficient Do{a) / Do{l) achieves a time-independent value [1^. Here, we compare the 
steady state values of DQ{a) / Do{l) obtained from Monte Carlo simulations with the theoretical predictions given by 
the first and second Sonine approximations. 

Let us consider first the dilute gas limit (0 = 0). Figure [1] shows the reduced diffusion coefficient Do{a) / Do{l) for 
mo/m = 1/8 and ao/a = 1/2 for disks {d = 2) and spheres {d = 3). According to the theory results obtained in 
the previous figures, one expects that in the Lorcntz gas limit (small values of the mass and size ratios) both Sonine 
approximations differ significantly for strong inelasticity (this means obviously that it is not possible that both of them 
simultaneously show good agreement with the exact solution of the problem). For the sake of completeness, we have 
also included the results recently derived for binary mixtures [25| from a new method based on a modified version of the 
first Sonine approximation which replaces the Maxwell-Boltzmann distribution weight function (used in the standard 
Sonine approximations) by the homogeneous cooling state distribution [25l . [isj . This new method partially eliminates 
the observed disagreement (for strong dissipation) between computer simulations [i^ and theoretical results for the 
heat flux transport coefficients [l^ . The explicit expression for the coefficient Dq obtained from the modified Sonine 
method is displayed in Appendix [Dl The simulation data corresponding to d = 3 for a > 0.5 were reported in Ref. 
[20| while those corresponding to d = 2 and d ~ 3 for a < 0.5 have been obtained in this work. It is quite apparent 
that while the second Sonine approximation agrees very well with simulation data, the standard and modiflcd first 
Sonine approximations fail for strong dissipation. This agreement is specially significant in the case of hard spheres. 
Thus, as expected, the Sonine polynomial expansion exhibits a poor convergence for sufficiently small values of the 
mass and size ratios. To assess the infiuence of density on these trends, the ratio Do{a)/Do{l) is plotted in Fig.[5]for 
mo/m = 1/5, ao/o- = 1/2 in the case of a moderately dense gas (</> = 0.2). As before, both first Sonine approximations 
underestimate the diffusion coefficient (the discrepancy being more important in the case of the standard than the 
modified first Sonine approximation) while the predictions of the second Sonine approach show an excellent agreement 
with simulation data in the whole range of values of the coefficient of restitution. All these results clearly confirm the 
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FIG. 7: (color online) Sketch of the segregation problem analyzed in Sec. IVII The small circles represent the particles of the 
dense gas while the large circles are the intruders. The BNE (RBNE) effect corresponds to the situation in which the intruder 
rises (falls) to the top (bottom) plate. 



accuracy of the second Sonine approximation, even for low values of a and small values of the mass and size ratios. 
Figure [S] shows that in the opposite limit of large values of the mass and size ratios (Rayleigh gas limit), the first 
and second Sonine approximations are practically indistinguishable for moderately large inelasticity (say for instance, 
a > 0.5) and both approaches provide a general good agreement with Monte Carlo simulations. However, at very 
high inelasticity, the second Sonine approximation very slightly underestimates the diffusion coefficient compared to 
the DSMC data and the first Sonine approximation. We have performed more series of simulations (not shown here) 
with different values of the ratios mo/m and cto/ct confirming similar trends as those in the figures shown in this 
Section for both cases toq/to < 1 and/or ao/a < 1 and mo/m > 1 and/or ao/a > 1 (the data are available to the 
reader upon request to the authors). 



VI. SEGREGATION BY THERMAL DIFFUSION 



As an application of the previous results, this Section is devoted to the study of segregation, driven by both gravity 
and temperature gradients, of an intruder in a granular dense gas. Segregation and mixing of dissimilar grains is one 
of the most interesting problems in granular mixtures not only from a fundamental point of view but also from a more 
practical point of view. This problem has spawn a number of important experimental, computational and theoretical 
works in the field of granular media [soj . Although several mechanisms have been proposed in the literature, the 
problem is not completely understood yet. Among the different mechanisms, thermal diffusion becomes the most 
relevant if the system resembles the conditions of a granular gas. In this case, kinetic theory tools have proven to be 
quite useful to analyze the motion of the intruder. A short previous analysis on this problem, but when the system 
is heated by a stochastic-driving force, has been reported in Ref. (3^ . 

Thermal diffusion is caused by the relative motion of the components of a mixture due to the presence of a 
temperature gradient. As a result of this motion, a steady state is finally reached in which the separating effect 
arising from thermal diffusion is balanced by the remixing effect of ordinary diffusion [5l[. Under these conditions, 
the thermal diffusion factor A characterizes the amount of segregation parallel to the temperature gradient. Our goal 
here is to determine A in a non-convecting (u = 0) steady state with gradients only in the vertical direction (z axis) 
for simplicity. In this case, A is defined as 

- A9^1nT = (9^1n(^) . (57) 

If we assume that gravity and thermal gradient point in parallel directions (i.e., the bottom plate is hotter than the 
top plate), then the intruder rises with respect to the ffuid particles if A > while the intruder falls with respect to 
the fluid particles if A < 0. If the impurity is heavier than the gas particles, the former situation is referred to as the 
Brazil-nut effect (BNE) while the latter is called the reverse Brazil-nut effect (RBNE) . A sketch of the geometry of 
the segregation problem studied here is given in Fig. [T] The key point here is that logically the segregation condition 
([57)) will depend on the mass flux transport coefficients of the intruder. We remind the reader of the features, some 
of them described in the Introduction, that the transport coefficients we use here are able to capture (compared to 
alternative coefficients in theoretical works on segregation by other authors [2^,[33|). These features are a consequence 
of the consistent development of the Chapman-Enskog theory for dense granular mixtures [l^ . Due to this, and in 
addition to the gain of accuracy with the use of the second order term in the Sonine expansion, we expect that our 
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FIG. 8: (color online) Plot of the thermal diffusion factor A versus the (reduced) gravity \g*\ for (f) = 0.2 and three values of 
the (common) coefficient of restitution a: a = 1, 0.8 and 0.5. 



expressions P5)) - (|T7)) will be useful in practical problems and applications of granular segregation [53|. We determine 
below the explicit expression of this segregation criterion. 

As said before, let us consider a steady base state with no flow (u = 0) and with a temperature gradient parallel to 
the direction of gravity (in this case, the z direction). According to Eq. (flBl) . the mass flux jo vanishes in this state 
(because u = 0) and there are no contributions to the pressure tensor except for those coming from the hydrostatic 
pressure term, i.e., Pij = pSij, where p is given by Eq. (|30p . As a consequence, the momentum balance equation ([9|) 
becomes 



dp 
dz 



n = -pg. 



Finally, the constitutive equation for the mass flux is given by Eq. ([T]) with V 
taking into account Eqs. ([T]) and ((57)) . the factor A can be written as 



dz- Using the fact that jo, 2 



(58) 
and 



A 



(3D'* -{p* +g*){Dl+D*) 
/3D* 



(59) 



where we have introduced the reduced transport coefficients D — {piy/noT)D , Dq — {viqU / pT)Dq, and D 
{moi'/noT)D. Moreover, p* = p/nT = 1 + 2'^-^{l + a)x''°U, 



(3^p* + 



and 



= l + 2'^-2(l + a)x(°V 



9' 



P9 



1 + 0— M^x^o); 



im 



< 



(60) 



(61) 



is a dimcnsionlcss parameter measuring the gravity relative to the thermal gradient. This parameter measures the 
competition between these two mechanisms (g and dzT) on segregation. 

As expected, the dependence of A on the parameter space of the problem is quite intricate. Regarding the depen- 
dence of thermal diffusion on gravity, we observe that for given values of the mass and size ratios, the coefficient of 
restitution and density, it is possible to switch between RBNE (A < 0) and BNE (A > 0) by changing the value of 

Pvity relative to the thermal gradient. This is a new interesting effect not captured in previous works on segregation 
, [lo, [13. As an illustration of this effect. Fig. [8] presents plots of the second Sonine approximation for A as a 
function of the (reduced) gravity for a three-dimensional system with </> = 0.2, ao/<7 = mo/m = 2 and three different 
values of the (common) coefficient of restitution a. It is apparent that, for the case analyzed here, the RBNE is 
dominant at small values of \g* \ while the opposite happens as the dimcnsionlcss gravity increases. 

The condition A = provides the segregation criterion for the transition BNE <^ RBNE. Since /3 and Dq are 
positive, then, according to ([55]) . sgn(A) = sgn{f3D^* — {p* + g*){DQ + D*)). As a consequence, the segregation 
criterion is 



PD^ 



{p* +g*){Dl+D*). 



(62) 
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FIG. 9: (color online) Phase diagram for BNE/RBNE for 4> = 0.1 in the absence of gravity and for two values of the (common) 
coefficient of restitution a. Points above the curve correspond to A > (BNE) while points below the curve correspond to 
A < (RBNE). The dashed line is the result obtained from the first Sonine approximation for a = 0.8. 
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FIG. 10: (color online) Phase diagram for BNE/RBNE for = 0.2 in the absence of thermal gradient {\g*\ — » oo) for three 
values of a: a = 1, 0.8 (dashed line) and 0.5 (dashed-dotted line). The dotted line refers to the results obtained by Jenkins 
and Yoon [2^ for an elastic system. 



As expected, when — m, ao ^ a and a = ao, the system (intruder plus gas) is monodisperse and the two species 
do not segregate. This is consistent with Eq. ([5^ since in this limit case D^* = + D* = so that A = for any 
value of a and (f). On the other hand, in the case of a dilute gas (0 = 0), one has [3 = p* — 1 and the condition ((62|) 
when g* ~ Q in the first Sonine approximation becomes 

^ = ^. (63) 

This segregation condition agrees with some recent results derived from the Boltzmann equation [3l|, [U ■ It must be 
remarked that, due to the lack of energy equipartition, the condition mo/m = Tq/T is rather complicated since it 
involves all the parameters of the system. 

We consider now dense systems. For the sake of concreteness, we assume that the intruder is larger than the gas 
particles (ctq > f)- Figure [S] shows a phase-diagram in the (mo/m, (To/o')-plane for a moderate dense gas {(j) = 0.1) in 
the absence of gravity and two values of the coefhcient of restitution. Also, for comparison the corresponding phase- 
diagram obtained from the first Sonine approximation is plotted for the case a = 0.8. We observe that, in the absence 
of gravity, the main effect of dissipation is to reduce the size of the BNE. This conclusion qualitatively agrees with 
the results derived in the driven gas case [1^. However, at a quantitative level, the influence of dissipation observed 
here is less important than the one obtained in the heated case. Moreover, although the first Sonine approximation 
reproduces qualitatively the trends of the phase-diagram, the former overestimates the predictions of the second 
Sonine approximation, especially when increasing the mass and size ratios. 

In some previous theoretical studies j29l . [sol . [52| , it has been assumed that the global temperature of the bed does 
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not vary with height so that the effect of dzT on segregation is neglected. In this limit {\g*\ oo), the condition 
(p^ reduces to Dq + D* = 0. The form of the phase diagram in the limit \g*\ — > oo is shown in Fig. [TO] for (j> ~ 0.2 
and three values of the coefficient of restitution. The result derived by Jenkins and Yoon [22] in the elastic case from 
a simple kinetic theory (namely, just a particular case from the perspective of the theory presented in this work) 
has also been included for comparison. In contrast to the case g* = 0, we observe that the RBNE regime appears 
essentially now for both large mass ratio and/or small size ratio. With respect to the influence of inelasticity, Fig. [TO] 
shows that the phase diagram is practically independent of the value of the coefflcient of restitution since the three 
curves collapse in a common curve. We also observe that our results differ from those obtained by Jenkins and Yoon 
[1^, especially for large size ratios. 



VII. CONCLUSIONS 



In this paper we have analyzed the mass transport of impurities in a moderately dense granular gas described by 
the inelastic Enskog kinetic equation. This is perhaps the simplest example of transport in a multicomponent granular 
gas since the tracer particles (impurities) are enslaved to the granular gas (solvent) and there are fewer parameters. 
Nevertheless, it involves, as we explained, many situations of practical interest in the study of granular gases. In 
the tracer limit, once the state of the solvent is well characterized, the mass flux jo associated with impurities is the 
relevant flux of the problem. To flrst order in the spatial gradients, jo is given by Eq. (|T]) where Dq is the kinetic 
diffusion coefficient, D is the mutual diffusion coefficient and is the thermal diffusion coefflcient. 

The main goal of this paper has been to determine these three transport coefficients as functions of the temperature, 
the density and the different mechanical parameters of the system, namely, the masses and particle diameters and 
the (constant) coefficients of restitution for the impurity-gas and gas-gas collisions. Like for elastic collisions li | , 
the coefficients Do, D and are given in terms of the solutions of a set of coupled linear integral equations [14|. 
A practical evaluation of the above diffusion coefficients is possible by using a Sonine polynomial expansion and 
approximate results are not limited to weak inelasticity. Here, Dq, D and have been determined in the first 
(one polynomial) and second (two polynomials) Sonine approximation and progress was possible here thanks to 
previous results obtained by using the (standard) first Sonine approximation [1^ for the full Navier-Stokes transport 
coefficients of polydisperse dense mixtures. The present study complements and extends previous works on diffusion 
in granular dilute 2^ and dense [s^ gases and provides explicit expressions for Dq, D and beyond the first Sonine 
approximation [15|. 

Comparison of the theoretical results derived for Dq, D and D^ between the first and second Sonine approximations 
shows significant discrepancies between both approaches for values of the mass ratio too/to and/or the size ratio oq/ a 
smaller than 1 while the quality of the first Sonine correction improves with increasing values of mo/m and cto/ct. 
These trends are quite similar to those previously found for tracer diffusion in an ordinary (elastic) dense gas (40| and 
in a granular dilute gas [2^. Moreover, to check to reliability of the different theoretical approaches, a comparison 
with Monte Carlo simulations of the Enskog equation for the coefficient Dq has been carried out for disks {d = 2) 
and spheres (d = 3). The comparison with simulation data shows the superiority of the second Sonine approximation 
over other approaches (the standard flrst Sonine approximation and a modifled version of the first Sonine correction 
recently proposed (25l.l48j). since the agreement of -Do [2] with numerical results is excellent, even for strong dissipation 
(see for instance Figs. H] and [5]) and very small values of the mass and size ratios. 

With respect to the segregation problem, which has of growing interest in the research community in the field, we 
have shown that the explicit knowledge of the three diffusion coefficients allows one to compute the thermal diffusion 
factor A. This quantity provides a convenient measure of the separation or segregation generated by a temperature 
gradient in a multicomponent system. According to the symmetry of the problem (sketched in Fig. [7]), when A > the 
intruder tends to climb to the top of the sample against gravity (Brazil-nut effect, BNE) while if A < the intruder 
tends to move at the bottom of the system (reverse Brazil-nut effect, RBNE). The understanding of the transition 
BNE/RBNE is of central interest in the field of granular matter mainly due to its practical/industrial importance. 
The analysis carried out here provides an extension of a previous analysis js^ for a heated dense gas in the flrst Sonine 
approximation. Our results show that the influence of dissipation on the phase diagram BNE/RBNE is much more 
signiflcant in the absence of gravity {\g*\ = 0) than in the opposite limit {\g*\ = oo). In fact, as Fig. 1101 shows, when 
the segregation of the intruder is essentially driven by gravity, the inelasticity of collisions has not discernible influence 
on the form of the phase diagram BNE/RBNE. From the discussion in previous Sections, we expect the seg regation 
criteria presented here to be more accurate in comparison with the criteria derived in other works p9l . [30l |. Future 
research work using a DSMC code adapted to the problem of segregation, molecular dynamics (MD) simulations, and 
eventually, experiments, will help us to test our theory (and previous alternative theories) in real problems. We are 
currently working on DSMC and MD simulations. 

An important issue is the usefulness of the expressions for the NS transport coefflcients derived here. As already 
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said in a previous work [25j , the NS hydrodynamic equations themselves may or may not be limited with respect to 
inelasticity, depending on the particular granular flow considered. In particular, tracer diffusion in the HCS at very 
low values of the coefficient of restitution is only possible for very small systems due to the spontaneous formation of 
velocity vortices and density clusters. Moreover, in most of problems of practical interest (such as steady states for a 
granular gas heated or sheared from the boundaries), the strength of spatial gradients is set by inelasticity so that the 
NS description only holds in the quasielastic limit [5a]. Nevertheless, in spite of the above cautions, the NS equations 
are still appropriate for a wide class of flows. Some of them correspond to the stability analysis of small perturbations 



of the HCS [57[ , supersonic flows past a wedge |58| and hydrodynamic profiles of systems vibrated vertically [59| where 
comparisons between theory and experiments have shown both qualitative and quantitative agreement for moderate 
values of dissipation (say for instance, a > 0.8). Consequently, the NS equations with the transport coefficients 
derived here can be considered still as an useful theory for a wide class of rapid granular flows, although more limited 
than for ordinary gases. 

One of the main limitations of the present study is its restriction to the tracer or intruder limit. This precludes the 
possibility of analyzing the influence of composition on the mass transport. The extension of the results derived here to 
finite mole fractions is an interesting open problem. Moreover, it would also be interesting to evaluate the expressions 
for the remaining transport coefficients of the mixture (shear viscosity, bulk viscosity, thermal conductivity, . . .) for 
a variety of mass and diameter ratios. This evaluation would allow one to assess the quality of the approximate 
Sonine method for solving the integral equations for the transport coefficients through a comparison with computer 
simulations. Previous results obtained for the shear viscosity coefficient jioj have shown a good agreement. Moreover, 
the knowledge of the full NS transport coefficients for a dense granular binary mixture allows us to determine the 
dispersion relations for the hydrodynamic equations linearized about the homogeneous cooling state. Some previous 
results [60j based on the Boltzmann kinetic equation have shown that the resulting equations exhibit a long wavelength 
instability for three of the modes. The objective now is to extend to higher densities this previous linear stability 
analysis for a dilute gas [60j and compare the theoretical predictions with MD simulations for the homogenous cooling 
state. We plan to carry out such studies in the near future. 
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APPENDIX A: FIRST ORDER VELOCITY DISTRIBUTION FUNCTION OF THE GAS 

The first order velocity distribution function /'^^' of the gas particles has the form Iv^ 

f^^^ = A-VT + C -Vn + T) ■.Vu + EV -u. (Al) 

Only the coefficients A and C are involved in the evaluation of the mass transport Jq^^ of the intruder. These quantities 
verify the linear integral equations 0, [l3| 



= A, (A2) 



2^ dV ' ' dn 
where J*^°) [X, Y] is the linearized collision operator 



1 ^(0) A • (VC) - n^A - ( J(°) [C, /(o)] + [/(°) ,e])^C, (A3) 



J(°) [vi|A,r] = x^"V'"' J dw^j dae{a ■ g)(S • g) [a-'X{V'l)Y{V'^) - XiY,)YiV^)] , (A4) 



and the inhomogeneous terms of the integral equations (jA2|) and (jA3P are defined by 

A (Y)--V—- (vf^°A --—f^°'> + -IC 



dY \ ■' 



(A5) 
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Q (V) = -V/ 
Here, the operator /CJX] is given by 



(0) _,-^-i_^mdp_( I d}nx^- 



f(0) 



(A6) 



The functions ^ and C are zero in the first Sonine approximation [Ta . [r7| . In the second Soninc approximation, 
these quantities are given by 



where 



and 



^(V) ^ -/M(V)a S(V), C(V) ^ -/m(V)c S(V), 



/M(V) = n(^) exp 



2T 



S(V)=('imF^-^T)V. 



(A7) 



(A8) 



(A9) 



Substitution of (|A7|) into Eqs. (|A2[) and (|A3|) gives a set of closed equations for a and c. Multiplication of these 
equations by S(V) and integration over V yields a set of algebraic equations whose solution is p^.fl7j 



(AlO) 



nTv 



aT^v 



2. AC 



AO) 



(All) 



where ly = na'^ ^ jm, C^* is given by 



and 



,(d-l)/2 



rf(d + 2) v^r(rf/2; 



X(")(l + a) 



^ + -(d + 8)(l-a) 



(A12) 



(A13) 



It must be remarked that all the above expressions have been obtained by neglecting some non-Gaussian contributions 
to the zeroth-order distribution /^°^. The influence of these non-Gaussian terms is only significant for quite extreme 
values of dissipation [49| . 



APPENDIX B: FIRST AND SECOND SONINE APPROXIMATIONS FOR THE MASS FLUX OF 

IMPURITIES 



In this Appendix we determine the transport coefficients Z?Oj ^ and associated with the mass flux in the first 
and second Sonine approximation. In this case, the functions Aq, Bo and Co are given by Eqs. respectively 
while A. and C are approximated by (|A7p . Let us start with the thermal diffusion coefficient , which is defined by 
Eq. (|20[) . To get it, wc substitute first Aq and A. by their Soninc approximations and (|A7|) . respectively, and 
then wc multiply the integral equation (j23p by rrio'V and integrate over velocity. After some algebra, the result is 

i^,-C^"))D^ + ^i.2ao^Z,, (Bl) 
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where 



Zi = i^sa 



1 - 



m mp 

and we have introduced the coUision frequencies 

_ 1 
dnoTo 



mop 



1 

2dp 



dvmoViJCo. 



d 



V/ 



(0) 



dvmoV- J^°^[/o,mV,/(o)], 



1^2 



(B2) 

(B3) 
(B4) 
(B5) 



The operator /Co,z[^] is defined in Eq. ((32l) . The colhsion integral appearing on the right-hand side of Eq. (jB2 
involving this operator has been evaluated in Rcf. [lH with the result 



dv \ ■' 



(B6) 



where to = ao/a and Mq = mo/(m + mo). If only the first Sonine correction is retained (which means ag = a = 0), 
the solution to Eq. (|B1|) is 



pmo 



mop 



2 m 



(B7) 



Here, ^'"^[l] denotes the first Sonine approximation to . To close the determination of up to the second Sonine 
approximation, we multiply now Eq. (|23|) by So(V) and integrate over velocity to get 



(.4-2C^"))ao + ^(.5-C("))i?^ 



where 



1 



Z2 = -avQ + — 



mo 



To d{d+2)noTl 
2 



dvS'o,i/Co,; 



(9V ^ 



f 6 = 



d[d + 


2) noT3 


2 


mo 


d{d + 


2) noT2 


2 


mo 




- 2) rioT^ 



-^|dvSo-Jr[/o,MSo,/(°)], 

L /■ rfvSo-J^l/o.MV, /(")], 
J 



;^/rfvso.jr[/rjMSo]. 



(B8) 

(B9) 
(BIO) 
(Bll) 
(B12) 



The collision integral of (jB9j involving the operator /Co,i is given by \\1 



^Jdv 5o,.(V)/Co,, 



— .(V/W) 



1 nMoT= 



.(1 + c.)^x[,°V(l + "0) [id + 2){Ml 1) 



+ (2d - 5 - 9ao)MoAf + (d - 1 + 3ao + 6a^)M2] + QM'^{1 + ao)^} 



(B13) 



where 7 = Tq/T is the temperature ratio and il/ = m/mQ + m. In reduced units and by using matrix notation, Eqs. 
(IBll) and (IB8I) can be rewritten as 



7^ '^2 

j/| - 2C* 



(B14) 
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Here, i^* = vi/v, D"^* = {Tni'/xQT)D^ , = T'^va^, a* = T'^va, and 

XI - - - 7) + \{l + u^fM,x^^U{^ + «o), (B15) 

X I ^ [(d + 2)(M2 - 1) + (2d - 5 - 9ao)AfoM 

+ - 1 + 3ao + 6ao)Af^] + 61/^(1 + ao)^} . (B16) 

The solution to Eq. (jB14p provides the explicit expression of the second Sonine approximation D"^[2] to D"^. It can 
be written in the form (j44p . where the dimensionless function H is 

The determination of the first and second Sonine approximations to the diffusion coefficients and D follows 
similar mathematical steps as those made before for . Here, only the final results will be provided. The kinetic 
diffusion coefficient Dq is obtained from the integral equation ([M)) by substitution of the Sonine approximation (|35p . 
In matrix notation, the coefficients _Do and 60 obey the equation 



where = {mlu/pT)Dn and = Tvb^. The solution to (|B18l) gives the second Sonine approximation to _Do- It 
can be written as Eq. where the dimensionless function F is 



F 



The corresponding matrix equation defining the coefficients D and cq is 

/ _ 1 ^* ^2^* \(D*\_(Y*- ^^vlc* 



(B19) 



where D* = (niQi' / noT) D , and Cq = uTvcq. Moreover, the inhomogeneous terms are given by 

X { [(d + %)Ml + {7 + 2d- 9ao)MaM + (2 + d + 3ag - 3ao)M^] 

+3M2(1 + aofe^ + [{d + 2)Ml + {2d - 5 - 9ao)MoM + (d - 1 + 3ao + 6al)M^] 9^ 

+ 2)0(1+0)}, (B22) 

where c* = nTvc, = mQT/niTQ is the mean-square velocity of the gas particles relative to that of the intruder 
particle and 

* ^ 7-^(C* - ^l)(^r - 7^^l«*) + (^1 - C){X^ - ^6*«*) .p.2o^ 

° i.2*(C-i^5*) + W-C)K*-2C) ■ ^ ' 

The second Sonine approximation D[2] can be easily obtained from Eq. (jB20[) and the result can be written in the 
form where G is 

^ ^ W ~ IC - 7^^2 (^2* - Kc*) {Y{ - 7^^|c*)"' ^^24) 
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Finally, in order to get the explicit dependence of the transport coefficients on dissipation, one still needs to compute 
the temperature ratio 7 = Tq/T. It is determined from the condition = C*j where Co is [Tsj 



Co = 



dT{i 



1/2 



(1 + ao) 



l-^(l + ^)(l + «o) 



(B25) 



APPENDIX C: COLLISION INTEGRALS 



In this Appendix we obtain the expressions for the collision frequencies v* . Except V2 and , the other quantities 
were already determined (isj for arbitrary composition. For the sake of completeness, we display now their explicit 
forms in the tracer limit (xi — *■ 0). They are given by 



dr(f) 



X^A'/a + ao) 



1 



1/2 



(CI) 



M-V)I2 



d{d- 



d-1 
(T \ (0) , 



3/2 



A-{d + 2)^B 



(C2) 



d(d + 2)r(f) \(7 



2^(d-i)/2 faY-' (0) , ^f0 



(C3) 



where 



TT 



(.-l)/2 ^-^-1 (o)M2 



3/2 r 



C+(d + 2): 



1 + e 



2ai- 



d + 3 



ao + d+l][d + 5+{d + 2)0] 



-D 



-M(l + 6) {X9-^[{d + 5) + (d + 2)61] [(11 + d)ao -5d-7] 
-0-^[2Q + d(15 - 7ao) + d^l - ao) - 28ao] - {d + 2)^(1 - ao)} 
+3{d + 3)X'^0-^[d + 5+{d + 2)0] + 2X0-^ [24 + lid + d^ + {d + 2)^6] 
+ {d + 2)e-^[d + 3 + {d + 8)6] - (d + 2)(1 + 0)e-^[d + 3+{d + 2)0], 



(C4) 



(C5) 



B = (d + 2)(l + 2A) + i\/(l + 6i){(d + 2)(l-ao)-[(ll + rf)ao-5d-7]Ae'"i} 
+3{d + 3)X'e~^ + 2M2 (^2a'a - ^-^"12 +d+l^0-\l + 0)^ 
-{d + 2)0-\l + 0), 



(C6) 



C = 2M2(1 + 0)^ (^2al - ^^Q'o + d + 1^ [d + 2 + (d + 5)61] 

-M(l + 61) {A[d + 2 + (d + 5)6l][(ll + d)ao - 5d - 7] 
+6'[20 + d(15 - 7ao) + d2(l - ao) - 28ao] + {d + 2)^(1 - ao)} 
+3(d + 3)A2[d + 2 + (d + 5)6*] - 2A[(d + 2)^ + (24 + lid + d^)e] 
+ (d + 2)6'[d + 8 + (d + 3)61] - (d + 2)(1 + 6i)[d + 2 + (d + 3)6*], 



(C7) 
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D = {d + 2)(2A ~e) + M(l + e) {{d + 2)(1 - ao) + [(11 + d)ao ~M- 7]A} 

d + 3 



-3((i + 3)A'' - 2M^ 2a, 



ai2 + d + 1 (1 + e)-^ + (d + 2)(1 + 61), 



(C8) 



Here, A = Mo(l — 7 ^). It must be noticed that Eqs. (jCip - (|C8p have been obtained by taking MaxweUian distributions 

for the reference homogeneous coohng state distributions Z*^"' and 

It only remains to evaluate the coUision integrals defining the collision frequencies and v^. To compute them, 
we use the property 



/ 



dvi/i(vi)J^") [Vi|F,G] = X^o^^~' I dvi I rfv2 F(vi)G(v2) 

X / dae{a-g){a-s) \h{v[)-h{Vi] 



with 

Use of this property in Eq. (jB4p gives 



V'l = Vi -i\/(l + ao)(?-g)^ 



dV 



siy(^) " xrM(i + a„)^?^+^ / dc. 



dc 



2 xe 



, d + 2 



(x • Ci), 



(C9) 



(CIO) 



(Cll) 



where c, = Vi/wo, x = g/vo, and wc have taken the MaxweUian approximation (|A8[1 for The integrals appearing 
in (jClOP can be evaluated by the change of variables {ci,C2} {x,y}, where y = Oci + C2, the Jacobian being 
(1 + 6')^'*. With this change the integrals can be easily performed and the final result is 



-1/2 



Similarly, the reduced collision frequency is given by 



dT(i) V" 



Mo 



(C12) 



(C13) 



APPENDIX D: MODIFIED SONINE APPROXIMATION 



The expression for the kinetic diffusion coefficient Dq derived from a modified version of the first Sonine approxi- 
mation recently proposed psi . l48j is displayed in this Appendix. This coefficient is given by 



7 



1 z^* 



Here, 



dr{d/2) 



X(°)(l-a2)(l + 



32 



(Dl) 



(D2) 



where e is the fourth cumulant of the gas distribution function f^^\ It measures the departure of from its 
MaxweUian form and its expression is [6ll | 



e{a) 



32(1 - a)(l - 2^2) 



9 + 24d - (41 - d)a + 30a2(l - a) ' 
The (reduced) collision frequency is given by [25l | 



dr\ 



2^(<i-i)/2 ^^Y- (0),,..^ /I 



1/2 



1 (3 + 46l)eo - 6*26 



16 



(1 + 0)^ 



(D3) 



(D4) 
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where gq is the corresponding fourth cumulant for the distribution function /q*^'' of the intruder. Its explicit form can 
be found in the Appendix of Ref . [B^I • 
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